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Abstract A triangulated spherical surface model is numerically studied, 
and it is shown that the model undergoes phase transitions between the 
smooth phase and the collapsed phase. The model is defined by using a 
director field, which is assumed to have an interaction with a normal of the 
surface. The interaction between the directors and the surface maintains the 
surface shape. The director field is not defined within the two-dimensional 
differential geometry, and this is in sharp contrast to the conventional sur- 
face models, where the surface shape is maintained only by the curvature 
energies. We also show that the interaction makes the Nambu-Goto model 
well-defined, where the bond potential is given by the area of triangles; the 
Nambu-Goto model is well-known as an ill-defined one even when the con- 
ventional two-dimensional bending energy is included in the Hamiltonian. 

Keywords Triangulated surfaces • Collapsing transition • Surface 
fluctuation • First-order transition 



1 Introduction 



Shape of membranes sensitively changes depending on certain specific envi- 
ronmental conditions such as flow fields, gravity, and thermal fluctuations 
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PP. Transformations of the shape, as well as the surface fluctuations, are 
typical to soft materials such as biological membranes [2]. Considerable 
efforts have been given so far to understand these phenomena through 
statistical mechanical treatments especially from the view point of phase 
transitions [HP]- 

Surface models for such phenomena are conventionally defined by using 
curvature Hamiltonians 0E1[7]. A well-known two-dimensional curvature 
energy is the so-called Helfrich and Polyakov Hamiltonian, which is rigor- 
ously defined by using the notions of two-dimensional differential geometry 
and plays a role for maintaining the surface shape [8"Il9"IfTU] . A linear bend- 
ing energy in the compartmentalized surfaces can also maintain the surface 
shape [TT | fT2 l fT3 p i4 p i5j. A unit tangential vector along the compartment 
defines the linear bending energy in those compartmentalized models [TT1 
[T^l[TnimT5] , while a unit normal vector of the triangles defines the two- 
dimensional bending energy of Helfrich and Polyakov in the conventional 
models |16 [ IT7 l lT8]. Thus, we see that the mechanical strength of the surface 
is always provided by objects defined within the surface geometry such as 
the linear bending energy and the two-dimensional bending energy. 

On the other hand, director fields are crucial to understand phenomena 
such as the main transition in liquid crystals including Langmuir monolayer 
[2] . The directors represent lipid molecules and can simply be described by 
the three-dimensional vectors on the surface. The chirality of membranes is 
also connected with the tilt of directors [19,20, 21 ,22, 23, 24 . The directors 
align to each others and become ordered at low temperature, while they 
become disordered at high temperature. It should also be noted that the 
directors, unlike the normal vectors of the surface, cannot be defined within 
the surface geometry. 

However, it is unclear at present whether the directors maintain the 
surface shape of membranes and what is the role of the director if it could 
maintain the surface shape. Therefore, interactions between the director 
fields and the surface are very interesting and still remain to be studied. 
We know that self-avoiding interactions can also maintain the surface shape 
against the collapse [25,26,27,28], however, the collapsing transition of self- 
avoiding surfaces is very time consuming for numerical studies. 

In this paper, we investigate whether the surface shape can be main- 
tained only by interactions between the directors and the surface. The 
problem we are interested in is whether or not the directors can provide 
mechanical strength to the surface for maintaining the shape. Two types of 
bond potentials are examined; the first is the conventional Gaussian bond 
potential and the second is the Nambu-Goto potential. It is well known 
that the surface model with the Nambu-Goto potential is ill-defined [29] . 
and the model is also ill-defined even when the curvature energy is included 
in the Hamiltonian. Therefore, we check whether the model with director 
is well-defined when the Nambu-Goto energy is assumed as the bond po- 
tential. Moreover, our interest focuses on the phase structure of the model 
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if the surface shape is geometrically well-defined. In that case, it is also in- 
teresting to see whether or not the phase structure of the model is different 
from those of conventional surface models defined by the above mentioned 
curvature Hamiltonians. 



2 Models 

The triangulated sphere, where the models are defined, is constructed from 
the icosahedron and is identical with those in [17] . By partitioning the 
icosahedron such that a bond of the icosahedron is split into i pieces, we 
have a triangulated surface of size N = 10£ 2 +2, in which 12 vertices are of 
coordination number q~5, and the remaining N — 12 vertices are of q = 6. 

The models are defined statistical mechanically and hence by the fol- 
lowing partition function: 

f/ n 

Z = J2 YldX ie xp[-S(X 7 d)}, S{X,d) = Si +bS 2 . (1) 

d J i=l 

The parameter b is the microscopic bending rigidity and is of unit kT, where 
k and T are the Boltzmann constant and the temperature, respectively. 
S(X, d) is the Hamiltonian, which is dependent on the variables X and d; 
X represents the three dimensional position of vertices and d represents 
a three dimensional unit vector, which will be defined below. The symbol 
/'n*Li dXi denotes that the center of mass of the surface is fixed in the 
integrations. Si and S2 are defined as follows: 

Si = J2( X i ~ Xjf, Sa^X^l-di-n^], (model 1), (2) 
and 

Si=J2 A *, & = X)2[ 1 - d «' n J(i)]. (model 2). (3) 

A i 

The bond potential Si of model 1 in Eq. Q is the Gaussian bond potential, 
which is defined by the sum of bond length squares, while Si of model 2 
in Eq.Q is called the Nambu-Goto energy, which is defined by the sum of 
the area of triangles A . 

The symbol di in 52 of Eqs.(|2]) and ([3]) is a three dimensional unit 
vector defined at the vertex i. We call dj as the director field or simply as 
the director. The symbols TO-jU) in 52 is a unit normal vector of the triangle 
j(i) surrounding the vertex i. The definition of 52 in model 2 is identical 
to that in model 1 . Since 52 is very similar to the bending energy of the 
first model in [15], we call 52 as the bending energy. The director field d^ 
and the unit normal vectors xiju) are shown in FigJT] 
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The difference between modei 1 and modei 2 is seen only in the defi- 
nitions of Si. We should note that the Nambu-Goto energy Si of Eq.© 
makes the model ill-defined if S2 is given by the conventional curvature 
Hamiltonian such that S2 = Yl(ij) (1 — n i ' n j)- ^ should also be noted that 
the director energy S2 in Eqs.f2]) and ([3|) is not defined within the surface 
geometry because the director is an external variable of the surface; it 
is not constructed without points outside the surface. For this reason S2 
is completely different from the bending energy of the first model in [18j , 
although the expressions of S2 in Eq. ((2]) and the bending energy in [18] are 
very similar to each other as mentioned above. 

We note that the models are symmetric under the transformations such 
that n — » — n and d — » d, where the normal vector n is chosen to have 
only one of the two-orientations of the surface. This symmetry implies 
the existence of the phase transition between two potential minima in the 
smooth phase. Although the orientation can not change from one to the 
other in the case of self-avoiding closed surfaces, it can change in our mod- 
els, which are phantom. This symmetric property is identical to that of 
the conventional curvature surface model on the phantom spherical sur- 
face, where the model is symmetric under n — ► — n. Therefore, the phase 
structure of the models in this paper is expected to have the same phase 
structure as the conventional curvature model. 

It should also be noted that the bending energy in this paper is different 
from the so-called elastic energy (l/2)fc t (nAd) 2 of Hclfrich in 5 , where n is 
the normal of the surface and d is the average orientation of the molecules, 
and kt is an elastic modulus. Since (nAd) 2 can also be written as sin 2 8 by 
using the angle 6 between n and d, then (n A d) 2 appears to be equal to 
1 — d • n in S2 of Eqs.([2]) and §5§ at sufficiently small 8. However, 8 is not 
always constrained to be sufficiently small in this paper. In fact, the angle 
in this paper has values in < 8 < it while d in [5] is allowed to have 
values only in < 8 < n/2; this is because the variable d in [5] is assumed 
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to be on the unit half-sphere while d in this paper is assumed to be on the 
whole unit sphere. 



3 Monte Carlo technique 

The dynamical variables X and d of the models are integrated out and 
summed over in the partition function of Eq.|T]). The integrations and 
the summations can be performed by the canonical Metropolis Monte 
Carlo (MC) technique on the triangulated surfaces. The random three- 
dimensional shift of vertices X — > X' = X+5X is accepted with the prob- 
ability Min[l, exp(— 5S)}, where SS = S 1 (new) — 5* (old). The random vector 
SX is chosen in a small sphere, whose radius is fixed at the beginning of 
the simulations so as to have 50% acceptance rate. The variable d can also 
be updated in almost the same technique as that of X. The new position 
d' is chosen on the unit sphere maintaining about 50% acceptance rate. 
N consecutive updates of X and those of d make one Monte Carlo sweep 
(MCS). 

The phase transitions are relatively strong rather than those in the 
conventional curvature models. For this reason, we concentrate on relatively 
small lattices of size up to N = 3612 in model 1 and N — 5762 in model 
2. The transitions are hardly computed on large sized lattices, because the 
correlation time becomes longer and longer with increasing N. 

The total number of MCS after sufficiently large thermalization MCS 
is about 1 x 10 8 ~4x 10 8 for the model 1 surfaces of size N = 812 and 
A/=1212, 8 x 10 s - 15 x 10 8 for those of iV=1962, 7V = 2562, and 7V = 3612. 
About 5 x 10 8 MCS is performed for the model 2 surfaces of size TV = 1692, 
and 10 x 10 s - 16 x 10 8 MCS for those of N = 2562, TV = 3612, and N = 5762. 



4 Results 

4.1 Snapshots, mean square size, and bond potentials 

Firstly, we show snapshots of surfaces and their surface sections of model 1 
and model 2 in Figs. [2a) -[2d). The surfaces are obtained in the smooth 
phase at the transition points b = 0.787 and b = 0.8545 of model 1 and 
model 2, respectively. Small pins or burs seen on the surfaces represent the 
director fields d^, which are of unit length. The directors are hardly seen 
on the surface of model 2, because the surface size of model 2 in Fig. [2jb) 
is quite larger than that of model 1 in Fig. [2(a). From the surface section 
of Fig. [2(c), we see that almost all directors turn inside the surface. This 
indicates that the orientation of the surface in Fig. [2a) is opposite to the 
initial one. The surface size can be characterized by the mean square size 
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(a) (b) 




m (A) 



Fig. 2 (Color online) Snapshots of the surfaces of size iV = 3612 of (a) model 1 
at 6 = 0.787 and (b) model 2 at b — 0.8545, (c) the surface section of (a), and (d) 
the surface section of (b). The mean square size is X 2 ~44 in (a) and X 2 ~331 
in (b) . Pins or burs sticking out of the surfaces are the directors d; . 

X 2 denned by 

i i 

where the symbol X in Eq.(|3]) is the center of mass of the surface. In fact, 
we have X 2 ~44 in Fig.^a) and X 2 ~331 in Fig.^b). 

The surface shape can also be reflected in the mean square size X 2 . 
Figures a) and[3Jd) show X 2 versus b of model 1 and model 2, respec- 
tively. The solid curves in Figs. Efa), OJb), EJJd), and[3{e) are drawn by the 
multihistogram reweighting technique [30]. We see that X 2 grows larger 
and larger against b with increasing N in both models. This indicates a 
collapsing transition between the smooth phase and the collapsed phase, 
although no discontinuity is seen in X 2 . To see the order of the transition, 
we plot the variance Cx 2 of X 2 in Figs.[3jb) and^e), where CV 2 is defined 

by 

C X , = (1/N)((X 2 -(X 2 )) 2 ). (5) 

The peak values C™2 X are plotted against N in Figs. |3fc) and [3^f) in a 
log-log scale. The straight lines in Figs.[3^c) and[3ff) are drawn by fitting 
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Fig. 3 (a) The mean square size X 2 vs. b of model 1, (b) the variance C X 2 vs. 
b of model 1, and (c) a log-log plot of the peak values C™2 X vs. N of model 1. 
(d), (e), (f) are those of model 2 corresponding to (a), (b), (c). The error bars 
in (a) and (d) are the standard errors, and those in (b), (c), (e), and (f) are the 
statistical errors. The solid curves are drawn by the multihistogram reweighting 
technique. 

the data to 

C™f x ~ N a , (6) 
where a is a critical exponent of the transition. Thus, we have 

(7 = 1.29 ±0.13, (model 1), 

(7 = 1.19 ±0.19, (model 2). (7) 

The value of a of model 1 is slightly larger than a — 1 and the one of model 
2 is almost identical to a — 1. Both results indicate that Cx 2 — *■ oo in the 
limit of TV — > oo, and therefore the order of the collapsing transition is 
considered to be first order from the finite-size scaling theory [3"l"ll3"2"ll3"3"] . 

The first-order nature of the transition can be seen more convincingly in 
the variation of X 2 against MCS at the transition point. Figures HJa) and 
Hfb) show X 2 vs. MCS, which were respectively obtained at the transition 
point b = 0.787 on the N = 3612 surface of model 1 and at b = 0.8545 on 
the A^ = 5762 surface of model 2. Two distinct states are clearly seen in the 
variations of X 2 ; one is the smooth state and the other is the collapsed state. 
The first-order collapsing transition can also be confirmed from a double 
peak structure in the distribution (or histogram) h(X 2 ) of the variation of 
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Fig. 4 (a) The variation of X 2 against MCS obtained at the transition point 
6 = 0.787 on the surface of size N = 3612 of model 1, (b) the variation at & = 0.8545 
on the surface of size TV = 5762 of model 2, and the normalized distribution (or 
histogram) h(X 2 ) of the variation X 2 of (c) model 1 and (d) model 2. 



X 2 . FiguresQJc) and[Hd) shows h(X 2 ) of model 1 and model 2, where the 
existence of the double peaks is clearly seen. 

We know that the curvature surface models undergo a first-order col- 
lapsing transition [171118] . therefore, the order of the transition remains un- 
changed if the curvature energy is replaced by the bending energy in this 
paper. However, the transitions in model 1 and model 2 seem rather strong 
than those of the curvature energy models. In fact, X 2 in the collapsed 
phase in Figs.[3{a) and[3Jd) is almost independent of N, and therefore, the 
Hausdorff dimension of the surface in the collapsed phase is expected to be 
H > 3, which is typical of strong transitions seen in phantom surface mod- 
els such as the tensionless model [T5] and the intrinsic curvature models 
[3"5ll36| . Figures a) and[5Jb) show X 2 vs. b of model 1 and model 2 on the 
N = 10242 and N = 16812 surfaces, which are relatively larger than those 
in Figs.^a) and^d). The total number of MCS is about 2 x 10 8 - 3 x 10 8 
in the smooth phase and 1 x 10 8 in the collapsed phase in both models, and 
these numbers seems insufficient for such large sized surfaces. However, we 
understand from Figs. [3a) and[Hb) that the transition is very strong. 
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solid lines are drawn to guide the eyes. 
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Fig. 6 (a) The Gaussian bond potential Si/N vs. b of model 1 and (b) the 
Nambu-Goto bond potential Si/N vs. b of model 2. The expected relation Si/N ~ 
3/2 is satisfied in both models. 



Finally in this subsection, the Gaussian bond potential and the Nambu- 
Goto bond potential Si/N are shown in Figs. EJa) and[Sfb), respectively, 
and we find that the simulations are successfully performed. From the scale 
invariance of the partition function, we have Si/N — 3(N — l)/(27V)~3/2 
if the center of mass of the surface is fixed. We find from Figs. EJa) and 
|6^b) that the expected relation S\/N ~3/2 is satisfied. 



4.2 Bending energy and surface fluctuations 

FiguresEta) and[7Jd) show the bending energy S2/(2Nb) versus b obtained 
in model 1 and model 2. The summation £\ m ^ ne definition of S 2 

in Eq. © or in Eq. © gives Y,iHjU) 1 = 2N b, where N B is the total 
number of bonds. This is the reason why 5*2 is divided by 2Nb in Figs. 
[T^a) andJTJd). Discontinuous changes of S 2 /(2N B ) are not so apparent in 
the figures just like X 2 in Figs.[3ja) and[3jd). The variance of S 2 is given by 
C S2 = (1/N){ (S 2 -(S 2 )f) and is shown in Figs. Efb) and^e). Peaks are 
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Fig. 7 (a) The bending energy S2/(2Nb) vs. b of model 1, (b) the variance Cs 2 
vs. b of model 1, and (c) a log-log plot of the peak values Cgl J ax vs. N of model 
1. (d), (e), (f) are those of model 2 corresponding to (a), (b), (c). The error bars 
in (a) and (d) are the standard errors, and those in (b), (c), (e), and (f) are the 
statistical errors. The solid curves are drawn by the multihistogram reweighting 
technique. 

seen in Cs 2 and indicate that the models undergo the transition of surface 
fluctuations. We show the peak values C!? ax vs. N in Figs.[7Jc) and[7Jf) in 
a log-log scale. The scaling behavior is observed such that C™ K ~ in 
both models, and we have 

H = 1.09 ± 0.20, (model 1), 

H = 0.84 ±0.25, (model 2), (8) 

where the fitting was done by using the largest four data in Fig. [7fc) and 
the largest three data in Fig. EJf). Thus, we see that fi ~ 1 in model 1, 
and therefore the order of the transition of surface fluctuation is considered 
to be of first order. In the case of model 2, fi is considered to be /i ~ 1 
within the error, and therefore the result is consistent with the first-order 
transition. We should note that the transition point b c (N), where Cg 2 has 
the peak, is almost identical to that for the collapsing transition, which 
is seen in Fig. [3] Thus, the transition of surface fluctuations is considered 
to occur at the same transition point of the collapsing transition in each 
model. 

Figures [H a) and Eljb) show the variations of S2/(2Nb) against MCS 
at the transition points of model 1 and model 2. S^/^-ZVb) in Fig. jSJa) is 
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Fig. 8 (a) The variation of S2/(2Nb) against MCS obtained at the transition 
point b = 0.787 on the surface of size N = 3612 of model 1, (b) the variation 
at b = 0.8545 on the surface of size N = 5762 of model 2, and the normalized 
histogram h(S2) of the variation S , 2/(2A r s) of (c) model 1 and (d) model 2. 

obtained at 6 = 0.787 on the iV = 3612 surface, while S 2 /{2N B ) in Fig.^b) 
is at b = 0.8545 on the N = 5762 surface. The histograms ^(S^) of model 
1 and mode 2 are shown in Figs. [8jc) and[8jd), respectively. We clearly 
see a double peak structure in both h(S2)- Thus, we have confirmed more 
clearly from the variations of Si/fiNs) and the histogram h(S2) that the 
order of the transition is of first order. 

The surface fluctuations are assumed to be reflected in the bending 
energy S2/(2Nb)- However, S2 includes the director field d^, which is not 
directly connected to the surface fluctuations. Therefore, it is interesting to 
see the bending energy S3 defined by S3 = n i ' n j)i where and 

rij are unit normal vectors of the triangles i and j, which have a common 
bond. Although S3 is not included in the Hamiltonian, S3 is considered to 
reflect the surface fluctuations. 

In order to see the transition of surface fluctuations, we plot S^/iVe 
in Figs. GDJa) andOJd). The corresponding variance Cs 3 defined by Cs 3 — 
(1/N)( (S3-{S 3 )) 2 ) is also plotted in Figs. Mb) and^e) against b. The 
peak values C^ ax against N are shown in Figs.[9^c) and[9jf) in the log-log 
scale. We find that the behavior of S3 /Nb is almost identical to that of the 
bending energy S2/(2N B ) in Figs.[7]Ja) and[7^d), and that the shape of the 
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Fig. 9 (a) The bending energy S3/NB vs. b of model 1, (b) the variance Cs 3 vs. 
b of model 1, and (c) a log-log plot of the peak values Cgi, ax vs. N of model 1. 
(d), (e), (f) are those of model 2 corresponding to (a), (b), (c). The error bars 
in (a) and (d) are the standard errors, and those in (b), (c), (e), and (f) are the 
statistical errors. The solid curves are drawn by the multihistogram reweighting 
technique. 



curves of Cs 3 is almost identical to that of Cs 2 in Figs.[7Jb) and [Tie). The 
exponent v defined by C!? ax ~ N v is obtained by fitting the data in Figs. 
E^c) and EK, and we have v= 1.13(18) and = 0.83(27) for model 1 and 
model 2, respectively. Both of v are consistent with the values of /j, in Eq. 
§8§ as expected. 

The variation of S 3 /N B vs. MCS is shown in Figs. QI|a) and[TI|b). 
Both of which are obtained at the same transition points where the varia- 
tions of S 2 /{2N B ) in Figs.E|a) andl^b) are obtained. We see that S 3 /N B 
changes between two different states corresponding to the smooth phase 
and the collapsed phase in both models, and the behaviors of the variation 
of S 3 /N B are almost identical with those of S 2 /(2N B ) shown in Figs.[S{a) 
and[Hb). This confirms that the transition of surface fluctuations is of first- 
order. Figures fTOf c) and flW d) show the normalized histogram h(S 3 ) of the 
variation S 3 /N B of model 1 and model 2, and we see in h{Ss) a double 
peak structure, which corresponds to that in h(S2) in Figs. Etc) and^d) 
and is consistent with the first-order transition of surface fluctuations. 

The correlation energy of the directors is defined by S4 = Yl(ij)0-~ ' 
dj), where denotes the sum over bond (ij) connecting the vertices i 

and j, and thus Yluj) 1 = N B . We checked whether or not the transition of 
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Fig. 10 (a) The variation of Sz/Nb against MCS obtained at the transition 
point 6 = 0.787 on the surface of size N = 3612 of model 1, (b) the variation 
at b = 0.8545 on the surface of size N = 5762 of model 2, and the normalized 
histogram h(S^) of the variation S^/iVe of (c) model 1 and (d) model 2. 



surface fluctuations is reflected in S4, and we reconfirmed from S^/Nb and 
the variance Cg 4 that the transition of surface fluctuations is of first-order. 
Although no figure of these quantities is shown, the behaviors of these 
quantities are almost identical with those of S2 in Figs. [7] and [8] and those 
of S3 in Figs. 151 and [TU1 The directors become parallel to each other and 
normal to the surface if the surface is sufficiently smooth, and therefore 
di are expected to be strongly correlated to each other in the smooth 
phase. On the contrary, are naturally expected to be uncorrclatcd on 
the fluctuated surfaces. Therefore, the transition of surface fluctuations is 
also reflected in S^/Nb- 



4.3 Scaling at b > b c 

The transition point b c is obtained from b c (N) in the limit of N — >oo, where 
b c (N) is evaluated from the peaks of Cx 2 in Figs.[3]Jb) and[3{e) or of Cs 2 
in Figs.ElJb) andCZle). Figures fTTTa) andQjIb) show b c (N) vs. 1/N in the 
linear scale, where b c (N) is obtained from the peaks of Cx 2 in Figs. EJb) 
and[3]Je). We have seen no clear difference between the transition point of 
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1/N 1/N 

Fig. 11 The transition point b c (N) vs. 1/N of (a) model 1 and (b) model 2. 
b c {N) are obtained from the peaks of C x i in both models. The straight lines are 
drawn by fitting the data linearly with respect to 1/N '. 



the collapsing transition and that of the transition of surface fluctuations. 
For this reason, we show only b c (N) obtained from the peaks of Cx^- 
The straight lines are obtained by assuming that b c (N) is proportional to 
1/N; b c (N) = b c + a(l/N), where the parameter b c is the transition point 
in the thermodynamic limit of the models. Thus, we have b c — 0.786 for 
model 1 and b c = 0.8539 for model 2. We should note that the values of b c 
are inconsistent with the transition points expected from X 2 in Figs. [5^a) 
and [DJb) , and the reason of this seems that the total number of MCS in 
the smooth phase was insufficient for the simulations of those large sized 
surfaces. 
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Fig. 12 The peak values C™J X vs. b/b c -l of (a) model 1 and (b) model 2 in a 
log- log scale, where b c = 0.786 in (a) and b c — 0.8539 in (b). 



The peak values C^2 X shown in Figs.[3fb) and[3^e) are plotted against 
b/b c ~l in Figs. [T2T a) and [T27 b) in a log-log scale, where b c = 0.786 and 
b c = 0.8539 are assumed. The straight line is the fitted one of the data such 
that C x 1 * x ~(b/b c -l)- a with exponents a = 0.92 ±0.16 and a = 0.97 ±0.15. 
These represent a conventional scaling property of C™^ at b > b c and 
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indicate that both models undergo the first-order collapsing transition. The 
results are also consistent with the values in Eq. (|7|). 
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Fig. 13 (a) The peak values C!? ax vs. bjb c —\ of model 1 and (b) the peak values 
C!? ax vs. b/b c — l of model 1 in a log-log scale, where b c = 0.786 in (a) and (b). 



C™ 2 ax and C*™ ax of model 1 are also shown in Figs. and H5Tb) 

against b/b c — l. The straight lines are drawn by fitting the largest three 
data in both of the figures. Thus, we confirm that the conventional scaling 
properties on the variances of 52 and S3 are almost consistent with the 
fact that model 1 undergo a first-order transition of surface fluctuations. 
However, in the case of model 2 we can not always see the scaling property of 
C™ 2 ax and C™ 3 ax in contrast to the case of model 1 in Figs.UHJa) and[l"3lTj). 
The reason of this seems only due to the low statistics of the simulations 
for model 2; the peak positions b c (N) corresponding to CE? ax and C™ ax in 
Figs.fTJc) and[5£c) are not always consistent with the plots for this scaling. 



5 Summary and Conclusion 

To summarize, we have investigated a possible mechanism that the direc- 
tors maintain the surface shape of membranes by using a spherical phantom 
surface model, which has no curvature Hamiltonian. The directors in the 
model are analogues of lipid molecules or some external molecules in mem- 
branes, whose shape is considered to be crucially dependent on the three- 
dimensional structure of the molecules. We focused our attentions on an 
interaction between the directors and the surface, and the phase structure 
of the model was studied by using the canonical MC simulation technique. 
We should note that the Helfrich Hamiltonian includes the mean curvature 
squared term and the Gaussian curvature term. The bending energy in this 
paper, as well as the conventional bending energy of the type 1— n • n, cor- 
responds to the mean curvature squared term in the Helfrich Hamiltonian. 
The Gaussian curvature term is eliminated from the Hamiltonian of the 
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model in this paper. This term seems relevant to the shape of such sur- 
faces with holes, however, it can be neglected in the case of closed surfaces 
because of the Gauss-Bonnet theorem [37] . 

Two-types of bond potentials are assumed in the model of this paper: 
The first is the standard Gaussian bond potential, which is included in 
the Hamiltonian of model 1, and the second is the so-called Nambu-Goto 
energy, which is defined by the area of the triangles and included in the 
Hamiltonian of model 2. Both model 1 and model 2 have the bending en- 
ergy, which describes the interaction between the directors and the surface. 
The shape of surface is maintained by this interaction. 

We found that both models undergo a first-order collapsing transition 
between the smooth phase and the collapsed phase. Moreover, a first-order 
transition of surface fluctuations occurs in the models at the same transition 
point of the collapsing transition. 

One remarkable result is that the model is well-defined even when the 
Hamiltonian includes the Nambu-Goto energy as the bond potential. This 
is confirmed from the numerical results of model 2. In the case when the 
Nambu-Goto energy is included in the Hamiltonian as the bond potential, 
the conventional curvature surface model becomes ill-defined [2T)j . 

The phase transitions seen in the models in this paper are relatively 
strong compared to those of the conventional models, although the order 
of the transitions is of first-order and hence identical to those of the con- 
ventional models. In fact, the collapsed surface at the transition point is 
completely collapsed in the models, while the collapsed surface of the con- 
ventional model is relatively swollen at the transition point. Consequently, 
the Hausdorff dimension H at the collapsed surface is greater than the 
physical bound, i.e. H > 3, in the models of this paper. This is in sharp 
contrast with the fact H < 3 in the collapsed phase at the transition point 
of the conventional curvature surface model |17j . 

We comment on the difference between the bending energy 1 — d ■ n in 
this paper and the elastic energy (n A d) 2 in [5]. The variable d in [5] has 
values on the unit half-sphere while d in this paper has values on the whole 
unit sphere. Therefore, the energy 1 — d ■ n in this paper is different from 
(n A d) 2 in [5], although both energies are almost equal to each other on 
sufficiently smooth surfaces. 

From the numerical results obtained in this paper, we conclude that the 
surface shape of membranes can be maintained by non-surface geometric 
object such as the directors, which interact with the surface. Important 
point to note is that the directors are the external variables of the surface. 
This implies that the directors are not always identified with the lipid 
molecules and allows us to speculate as follows: If some external objects 
could be embedded in the membrane so as to have the assumed interaction 
with the membrane constituents, the surface shape can be controlled. 

It is interesting to study the phase structure of the fluid surface model 
with the directors. Correlations between the directors can be assumed as an 
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energy term for maintaining the surface shape; we expect that the surface 
shape is maintained by the correlation energy. 

This work is supported in part by a Grant-in- Aid for Scientific Research 
from Japan Society for the Promotion of Science. 
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